Vibrational Sidebands and Kondo-effect in Molecular Transistors 
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Electron transport through molecular quantum dots coupled to a single vibrational mode is stud- 
ied in the Kondo regime. We apply a generalized Schrieffer- Wolff transformation to determine 
the effective low-energy spin-spin- vibron-interaction. From this model we calculate the nonlinear 
conductance and find Kondo sidebands located at bias-voltages equal to multiples of the vibron 
frequency. Due to selection rules, the side-peaks are found to have strong gate-voltage dependences, 
which can be tested experimentally. In the limit of weak electron-vibron coupling, we employ a 
perturbative renormalization group scheme to calculate analytically the nonlinear conductance. 

PACS numbers: 73.63.Kv, 73.23.Hk, 72.10.Fk, 72.15.Qm, 05.10.Cc 



In recent years, the study of transport in meso- 
scopic systems has branched into investigations of single- 
electron devices based on single-molecule transistors 0, 
1^1 IE ^ 111' Of particular interest is the possibility of 
combining electronics with mechanics, such that the vi- 
brational or configurational modes of the molecule are 
coupled to its charge state. A number of interesting 
issues have already been addressed in this new field of 
nano-electro-mechanics. Firstly, it was shown by Park et 
al. that quantum mechanical behavior of the center of 
mass oscillation of a Cgo can be excited by the tunneling 
electrons and a series of assisted steps were observed in 
the current. Similar structures have later been observed 
in a number of other experiments. These steps were asso- 
ciated with simple Franck-Condon physics, which means 
that the tunneling rates are modified by the overlaps of 
the initial and final states of the oscillator. 

It is well established H, 0] that single particle reso- 
nance tunneling is not destroyed by the electron-vibron 
coupling, but instead the resonance breaks up into a num- 
ber of vibron sidebands. The question remains, though, 
whether more intricate many-body effects, such as the 
Kondo resonance, also cooperate with the electron-vibron 
coupling to form "Kondo sidebands" . The usual Kondo 
resonance has been observed in several molecular de- 
vices H, at unusually high temperatures, and recent 
experiments ^ on Ceo, ^nd Co based transistors have 
revealed marked sideband resonances, which were sug- 
gested to arise from the interplay of a Kondo resonance 
with a vibrational mode. 

In this letter we demonstrate that, in contrast to se- 
quential tunneling, which is suppressed by the Franck- 
Condon overlap factors, the Kondo resonance remains 
intact well inside the Coulomb blockade valley. In fact, 
the electron-vibron coupling is predicted to enhance the 
exchange-coupling and thereby the Kondo-temperature. 
Maintaining the quantum coherence of vibrons, we show 
that the Kondo-resonance breaks up into a series of vi- 
bron sidebands. Moreover, we demonstrate that parity 
selection rules prohibit all sidebands at odd multiples of 



the oscillator frequency, when tuning the gate-voltage to 
the particle- hole symmetric point. 

Assuming the energy-level spacing on the molecule to 
be much larger than the charging- energy, the system may 
be described by the Anderson-Holstein Hamiltonian 

a,k,cT 

+ ^ {takdlcaWa +ii-c.) + XuJo{b + b^)nd. (1) 
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where 



ikcr 



and are creation operators for electrons 
in the left and right conduction bands {a = L, R) and 
on the molecular quantum dot, respectively, n^a — dl.da-, 
n-d = + Udi and S^ak = Ek — /^q • The vibrational mode 
of the molecule is created by the vibron operator b^ , and A 
denotes the dimensionless coupling strength. Describing 
the molecule as a quantum-dot, we have Sd = {l~2J\f)Ec 
and U — 2Ec, in terms of the charging energy Ec and 
the mean occupation number Af = CgVg/e, determined 
by the gate voltage, Vg, and the capacitance to the gate, 
Cg. Unless explicitly stated otherwise, we shall hence- 
forth work in units where e = h = ks = 

Following Lang and Firsov , the electron-vibron cou- 
pling in the Hamiltonian ^ is eliminated by the unitary 
transformation H' = e^^" He""^^" , with Sp = iX{b—b^)nd: 



H' = CakC^kcrCakfT + Ed^d + iOob^ b + U'ud^ridl 

+ $](tcke^^'''^'^4ccka + h.c.), (2) 
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where = — A^cjq and U' = U — 2X'^ujo. We now con- 
sider the weak-tunneling limit, r^k — 27riV(0)|tQ.kP <C 
nim{—e'd,e'd + U'), where iV(0) denotes the conduction 
electron (ce) density of states. In this limit, a generalized 
Schrieffer- Wolff transformation, devised by Schiittler and 
Fedro may be used to eliminate all first order terms 



2 



in ta. To this end, we introduce the generator 

a, k, (7,77 



(3) 



where = (1 — t])/2 + r/Udw, with r/ = ±1,(t = —a, 
and Caka,, - dte-*(^°''"-«+)*e-^(*), with ^„kr, = 
Cak - £rf " (1 + v)U'/2 and = A(e-"^°*6 - e'"^"*^^). 
Applying the transformation H" — e^^" H'e^'^^^' and ex- 
panding to second order in ta, one finds that H" = 
Hq + i?spin + -f^dir + -f^pair ■ We neglect the renormalization 
of the kinetic energy term in Hq and, restricting to the 
regime of single occupancy, i.e. A/" « 1 and Xloq <^ Ec 
{U' > 0), -ffpair vanishes. The potential scattering term, 
i/dir, is omitted since it leads to no logarithmic singular- 
ities, and altogether we obtain the effective Hamiltonian 



H" = ^ CakC^k(T"^ak<T + LUob^b 

(4) 



a,k,a" 



a,k,cr;a',k',fT' 



where S = ^d^/To-'crrfa denotes the local spin-1/2, and 

Ja,k;a',k' = C'k'*"k [(AT^k ~ ^ak) + i^a'k' ~ ^I'k'V] 

with = iJ^°^te"*('^°'"'"*+)*+'''^'' sin(i^ot)gA(o)-A(t)^ 

In this effective Kondo-model, the exchange-coupling JJ 
incorporates the dynamics of the vibron through the dis- 
placement operator e^. In the vibron number-state ba- 
sis it is therefore convenient to introduce Franck-Condon 
factors /„'„ = {n'\e~^^°'^\n) 0| , which allows us to write 
the matrix-elements of J in the more transparent form: 



'^a',k';a,k = |JlQ',k';a,k|'^ 
oo f 

^a'k'^CKk ^ ^ ^ Jinn' ft 



(5) 



m=0 



1 



fn'mfn 



Cak - £- + (to - n')ujo 

1 

£,a'-k' - e- + (to - n)ujQ 
1 

^Qk - e+ - (to - n)ujo 
1 



Ca'k' - e+ - (m - n')LLiQ 



valid for ^ak, Ca'k', Jt-wo, n'wo ^ min(e+,— e_), where 
e_ = and e+ = e'^ + U' are the energies of intermedi- 
ate, empty, or doubly occupied states. In this sum, the 
energies of intermediate vibron states |to) shift the energy 
denominators and the Franck-Condon factors determine 
the overlap between initial and final vibron states with 
intermediate states of the oscillator shifted by \/2Mo, 
where £o is the characteristic oscillator-length. 

Since q fn'mfnm — ^n'n and fn'n ^ ^n'n for 

A 0, the usual exchange-coupling, Ja>a — '^ta'ta/ Ec, 
is recovered in either of the limits ^ or A — s- 
0. More generally, J" " may be represented as an 



asymptotic power series as Ec/ojq oo, with lead- 
ing terms J^,'^ oc Ja'a(Awo/£^c)'"'~"'- In terms of the 
incomplete Gamma function, 7(0;, a;), one has = 
J„.„e-^'(£;cM)E„=±(-A')''^"/""7(»7e^/^o,-A2), or 
simply J"° « Ja'a[l + (Awo/Sc)^] for Xluq < Ec and 
A/" = 1, as found earlier in Ref. |9| . Staying well inside the 
Kondo-regime, any finite A thus leads to a slight enhance- 
ment of j"" , and thereby of the associated Kondo temper- 
ature, Tk ~ De-i/^^o)-^™ {2D being the ce bandwidth). 
This is essentially due to the vibron induced reduction of 
|£±|. Note also that the (cotunneling) amplitude 7°" in- 
volves only virtual shifts of the oscillator and therefore no 
Franck-Condon overlap reduces its magnitude. This is in 
sharp contrast to the resonant (sequential) tunneling am- 
plitude, relevant outside the Coulomb-blockade regime, 
which involves real excitations of the oscillator and is 
thereby reduced by a factor of |/ooP 'S, ^• 

We now consider the case of strongly asymmetric and 
momentum independent tunneling amplitudes, ^ 
r^. The current traversing the molecule from left to 
right is then given simply as llj 



I 



(6) 



From the equations of motion for the Hamiltonian l|T|l. 
the local density of states is found to be related to the 
ce T-matrix as Im[^;;^''^(a;)] = \tc,\~'^lT[v[T'Z^{uj)], and the 
latter can now be obtained using the effective Hamilto- 
nian To third order in Ja,a', we find that 



Im 



[T^;^(f^)] = -5,,,^7V(0)(l - e-o/^)(l + e-"/^) 



Imn 
aia2 



x9{D -\n + {n- l)ujo\) \ S„uSa,a, + A^(0) J^^ 



In 



D 



Q + (n — m)ujQ 



In 



D 



+ (m — l)ujQ 



,(7) 



with the shorthand notation In |_D/a;| = In \D/y/x'^ + T^j. 
In the asymmetric limit considered here, we take fj,L = 
and bias the right lead to fXR = —V, leaving the position 
of the molecular energy levels unaffected. For T ^ ojq <^ 
D, the differential conductance is then obtained from 
Eq. ® as G{V) = ~^{rn/TL)NiO)j:^lm[TllieV)]. 
From Eq. ||7J), the differential conductance appears to 
diverge as ln(D/r) at voltages corresponding to mul- 
tiples of the oscillator frequency, reflecting the onset 
of a Kondo-effect assisted by coherent vibron-exchange. 
In Fig. n the upper panel shows a gray-scale plot of 
d^I/dV^ as a function of bias- voltage V and mean oc- 
cupation number (gate- voltage) M — CgVg/e. The lower 
panel shows three cuts revealing the side-band resonances 
on the flanks of the central zero-bias resonance. 
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FIG. 1: Upper panel: d'^I/dV^ vs. bias and gate voltage, for 
= 3, N(Q)\tL\^ = O.Iljo, D ^ Ec ^ Sljq, and T = O.OIluo. 
Black/white indicates large negative/positive values. Lower 
panel: Conductance vs. bias voltage for three values of Vg 
corresponding to the vertical black lines (a,b,c) in the upper 
panel. The lower curve (a) corresponds to the ph-symmetric 
point TV = 1 - X'^ujo/Ec- 



By tuning the gate- voltage to TV = 1 — X'^ujq/Ec one 
reaches the particle-hole (ph) symmetric point where 
e_i_ = — and using the general symmetry = 
(-l)l"'-"l/„„- one finds from Eq. © that J"'" oc [1 + 
(-1)1"'-"!], implying that aU spin-exchange processes in- 
volving the emission or absorption of an odd number of vi- 
brons are prohibited at this particular gatc-voltage. This 
parity selection rule, reflecting the inversion symmetry 
of the Kondo-Hamiltonian at low energies, has im- 
portant experimental bearings, since it predicts that all 
side-band resonances in the differential conductance lo- 
cated at voltages equal to odd multiples of loq, must van- 
ish when tuning the gate- voltage to the symmetric point, 
corresponding to the midpoint of the Coulomb-blockade 
valley. This is apparent in Fig. ^ where the conduc- 
tance peak at V = loq disappears as Af approaches the 
symmetric point (cf. curve a). Note, however, that any 
appreciable vibron modulation of the tunneling ampli- 
tudes will break the inversion symmetry ,^10] and thereby 
destroy this selection rule. 

The logarithmic divergences appearing in third or- 
der perturbation theory call for a resummation of lead- 



ing logarithmic contributions to all orders. This is 
done using the perturbative renormalization group (RG) 
method for frequency dependent couplings developed in 
Ref. Parametrizing the dimensionless couplings, 

gn'n = N{Q)J2^, by the total energy of the ingoing 
conduction electron and vibron-state, we arrive at the 
infinite hierarchy of coupled (1-loop) RG-equations: 



d\nD 



[5«'m(0)gm„((m - n)wo)e 



LO-\-{n—rri)ujQ 



((n' — r7i)a;o)5mji(O)0w+(m-n')'-'o] i (8) 

with = Q{D — \uj\). We shall restrict our attention 
to the ph-symmetric point and assuming that Xluq <C Ec 
we may truncate this hierarchy and consider merely the 
lowest four coupled equations involving goo, go2 ~ .920 
and §22 ■ The solution to this reduced set of equations is 
characterized by the parameters 



922 — goo 

500522 - 920 ' 



a 



252 



.g00522 - 520 



(9) 



and Tk = De-2/(»"+9^^+V(s22-9oo)^-H4s2j^ ^^^^^ ^^^^ ^ 

9n'niD;2nujo). All three parameters are invariant un- 
der the perturbative RG-fiow from the initial cut-off 
Do down to 13 = 2(jJo- At scale Dq, we have a,S ~ 
{XuJo/Ec)^ /goo, and therefore our truncation of Eq. © 
remains valid throughout the RG-flow roughly when 
max(Q!, (5)/ ln(r/r*) <^ 1 (see below). Staying within 
the perturbative regime, we assume that wo 3> T ^ Tk- 
We first solve the RG-eqs. for the constant coefficients 
gn'n and the frequency dependent renormalized couplings 
are then obtained simply by integrating Eqs. (jHl): 



1 



500 



a 



1^=0,1 



ln(|^|/r*) 81n2(2^o/T*) 
1 



(10) 



1 



520 (w) = jJ2 



!^ = 0,l 



ln(|w - 2t/wo|/r*) ln(2wo/T*) 
1 

ln^(max(2a;o, 1^^ — 2i/[x)o|)/r*) 
1 1 



(11) 



ln(2wo/T*) \lii{\uj - 2iyujo\/T*) ln(2wo/T*) 



withT* = r^(r^/^o)(v'^W+5)/(21n(2.^o/T^) + x/^W-A-) 

and Qu{i^) — 6(2wo — |tf — 2voJo\)- Note that in Eqs. 
(|10lll|l we have retained only terms which will contribute 
to order max(a, (5)Vln^(r/T*) in the conductance. 

While the logarithmic singularities at a; = are cut 
off by temperature, those ai uj — 2u!o will instead be 
contained by a/T^ + 7^, with 7 given by the transi- 
tion rate from vibron-state |2) to |0). Using the golden 
rule with the renormalized coupling 1720 (w), we find 7 « 
711000^ / {Aln^ (2ll!o /Tk)) ■ Similarly, the broadening of the 
vibron states, induces a broadening of the step-fimctions 
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FIG. 2: Conductance versus bias voltage using perturbatively 
renormalized couplings. Gate voltage is tuned to the ph- 
symmetric point and a = 4.30, S = 7.54, Tk = 1.51 x lO^^'tJo, 
T* = 7.18 X 10"*a;o, and T = 0.05(Jo, i.e. a/ln(T/r*) = 0.32 
and 5/ln(T/T*) = 0.56, corresponding to bare parameters: 
= 4.5, N{0)\tL\'^ = O.Iluo and D = Ec = Swq. Insets show 
the renormalized couplings (;oo(<^) and g2o{ui), as well as a 
zoom in the conductance curve showing the satellite peak on 
a separate conductance scale but on same voltage scale. 

in Eqs. H1UI11|) by T or ^JT'^ + 7^ for steps near w = 
or a; = 2ijJo, respectively. 

The renormalized conductance is now obtained by in- 
serting in the formula G{V) = (2eV/i)(rfl/rL)(37rV4) 
[gmie-Vf + Y.u=± 'd{veV - 2u:Q)g2o{veVf], and indeed 
when expanding this result in bare couplings, we re- 
cover the result obtained by expanding Eq. iQ to order 
(Xluo/Ec)^- Including the broadening in both couplings 
and 8, the result is plotted in Fig. [3 This figure bears 
a certain resemblance to the curve obtained in Ref. [l^ 
at T <C Tpc, using an entirely different method. 

The possible decoherence effects which will arise when 
coupling the oscillator to phonons within the leads re- 
main an open question. It is straightforward to generalize 
the Schrieffer- Wolff transformation applied here to a sys- 
tem where the molecule-oscillator is coupled to a separate 
bath of oscillators. However, even determining the effects 
on the leading logarithms for a given Q-factor involves 
a rather involved cumulant expansion. We expect the 
Kondo-effect to be more pronounced when dealing solely 
with mtra-molecular vibrations, since these have been 
demonstrated to have particularly large Q-factors Q . 

In the case of nearly symmetric couplings, we can 
no longer assume the oscillator to be in equilibrium 
with the conduction electrons of one specific side of the 
i unction [ist. In Hne with the findings of Ref. 0,^^, we 
expect that nonequlibrium effects may in fact serve to 
enhance the Kondo side-peaks. 

In conclusion, we have demonstrated the viability of an 
inelastic Kondo-effect carried by coherent vibron-assisted 
exchange tunneling, which can be observed as Kondo- 
sidebands in the nonlinear conductance. In contrast to 
the case of an applied microwave field , the zero-bias 
resonance is not suppressed by the vibronic coupling, and 



it may therefore be difficult to discern the satellites from 
the background conductance. Nevertheless, even with 
very weak satellites, it should be possible to track their 
dependence on Vg (possibly in a plot of d'^I/dV'^) and 
thereby test our prediction that satellites at odd mul- 
tiples of 070 are strongly reduced at the ph-symmetric 
point. Faint sidebands to a zero-bias Kondo peak have 
indeed been observed in the recent experiment by Yu et 
al. 5j. In contrast to the findings reported here, how- 
ever, both satellite peaks and dips were observed in the 
nonlinear conductance. 
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